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Abstract 

High water use efficiency (WUE) can be achieved by coordination of biomass accumulation and water consumption. WUE is 
physiologically and genetically linked to carbon isotope discrimination (CID) in leaves of plants. A population of 148 
recombinant inbred lines (RILs) of sunflower derived from a cross between XRQ and PSC8 lines was studied to identify 
quantitative trait loci (QTL) controlling WUE and CID, and to compare QTL associated with these traits in different drought 
scenarios. We conducted greenhouse experiments in 2011 and 2012 by using 100 balances which provided a daily 
measurement of water transpired, and we determined WUE, CID, biomass and cumulative water transpired by plants. Wide 
phenotypic variability, significant genotypic effects, and significant negative correlations between WUE and CID were 
observed in both experiments. A total of nine QTL controlling WUE and eight controlling CID were identified across the two 
experiments. A QTL for phenotypic response controlling WUE and CID was also significantly identified. The QTL for WUE 
were specific to the drought scenarios, whereas the QTL for CID were independent of the drought scenarios and could be 
found in all the experiments. Our results showed that the stable genomic regions controlling CID were located on the 
linkage groups 06 and 13 (LG06 and LG13). Three QTL for CID were co-localized with the QTL for WUE, biomass and 
cumulative water transpired. We found that CID and WUE are highly correlated and have common genetic control. 
Interestingly, the genetic control of these traits showed an interaction with the environment (between the two drought 
scenarios and control conditions). Our results open a way for breeding higher WUE by using CID and marker-assisted 
approaches and therefore help to maintain the stability of sunflower crop production. 



Citation: Adiredjo AL, Navaud 0, Muhos S, Langlade NB, Lamaze T, et al. (201 4) Genetic Control of Water Use Efficiency and Leaf Carbon Isotope Discrimination in 
Sunflower (Helianthus annuus L.) Subjected to Two Drought Scenarios. PLoS ONE 9(7): e101218. doi:10.1371/journal.pone.0101218 

Editor: Carl J. Bernacchi, University of Illinois, United States of America 

Received January 16, 2014; Accepted June 4, 2014; Published July 3, 2014 

Copyright: © 2014 Adiredjo et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: The work was supported by a French Government Scholarship (Bourse du gouvernement Francais, BGF) and a co-funding by Directorate general of 
Higher Education, Ministry of Education and Culture, Republic of Indonesia (Beasiswa Luar Negeri, BLN). The funders had no role in study design, data collection 
and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* Email: grieu@ensat.fr 



Introduction 

Water use efficiency (WUE) as a breeding target can be defined 
as the ratio of biomass production to water consumption. Breeding 
for WUE and drought-resistant crop varieties has been a critical 
area of agricultural research worldwide [1-3]. Substantial efforts 
have been devoted to identifying and selecting for morphological 
and physiological traits that increase WUE and yield under rain- 
fed conditions [2,4-5]. In field conditions, water consumption is 
usually difficult to determine. Nevertheless, WUE can be 
represented by measuring leaf carbon isotope discrimination 
(CID) [6-7]. Because the CID has been demonstrated to be a 
simple but reliable measure of WUE, the negative correlation 
between them has been used as an indirect method in selection to 
improve WUE [8-10]. The principle mechanisms underlying the 
variation of CID act through variation in the intercellular C0 2 
concentration (c;) maintained in leaves [6]. The value of c, is 
determined through the coordinated regulation of carboxylation 



capacity (photosynthesis) and stomatal control of leaf diffusive 
conductance (transpiration regulation) [6-7]. 

Genetic variation underlying quantitative traits, such as WUE 
and CID, which are generally under considerable environmental 
influence, is governed by quantitative trait loci (QTL) [11-14]. 
QTL mapping provides a starting point in breeding programs 
[15-16] and for cloning of the causal mutations by fine mapping. 

QTL mapping of WUE is rarely reported. Four QTL associated 
with WUE have been identified in soybean [17]. The inheritance 
of WUE has been studied using simple sequence repeat (SSR) 
markers in alfalfa [18]. In contrast, QTL mapping of CID has 
been reported by numerous authors. The first QTL identified for 
CID was reported by Martin and Nienhuis [19]. These authors 
identified four QTL associated with CID in tomato. Since that 
time, QTL for CID have been identified across a wide range of 
species, for example in cotton [20], rice [21], barley [22], 
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Figure 1 . Principles of the water treatments used in this study. (A) In experiment 201 1 , three replicates (each of 1 50 plants) were subjected to 
progressive water-stress by water withholding from 1 to 31 DAE. In this experiment a control replicate (150 plants) was watered to maintain non- 
stressful conditions (SWC = 30%). (B) In Experiment 2012, two replicates (each of 150 plants) were maintained at in stressful conditions SWC= 16% 
from 1 to 23 DAE whereas two other replicates (each of 150 plants) were irrigated to maintain non-stressful conditions (SWC = 30%). DAE: day after 
emergence. 

doi:10.1371/journal.pone.0101218.g001 



Arabidopsis [23], and in wheat [24]. However, to our knowledge, 
QTL of WUE and CID in sunflower have never been reported. 

Most of the work identifying QTL of WUE and CID has been 
done in well-watered conditions, with only one study in a drought 
situation. There is no report on the QTL identification of WUE 
and CID of crops subjected to different scenarios of water deficit 
establishment. 

The objectives of the present study are to identify QTL 
controlling WUE and CID in a population of RILs of sunflower, 
and to compare QTL associated with these traits in a dual drought 
scenario: (i) a progressively water-stressed establishment and (ii) a 
stable water deficit treatment. We are interested in providing new 



insights into the genetic architecture of WUE and CID, and in 
contributing to the potential of sunflower breeding by improved 
WUE. 

Materials and Methods 

Plant materials 

A population of 150 recombinant inbred lines (RILs) was used 
in two experiments. A population of these RILs was named INEDI 
and was obtained by single seed descent (self-pollination to at least 
F8) from a cross between XRQand PSC8 [25]. 
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Experiments and trait measurements 

Two experiments were conducted in spring 2011 (Exp. 2011) 
and in spring 2012 (Exp. 2012) under quite similar weather 
conditions. Plants were grown in a greenhouse at the INRA 
Auzeville station, Toulouse, France (43°31 '46,94" N; 1°29'59,71" 
E). Greenhouse air temperature was set at 25/18±2°C (day/night) 
and relative humidity was about 55-75±5%. 

Three seeds per genotype were sown in a pot (volume: 2 liters) 
at the beginning of the experiments. The pots contained a mixture 
of 50% soil (collected from the field), 30% organic matter and 20% 
sand. These pots were arranged on 100 balances (maximum 
capacity 30 kg, precision 2 g, model SXS, GRAM, Spain), with six 
pots per balance (total pot number in greenhouse was 600). Each 
pot was then covered with a 3 mm layer of polystyrene sheet with 
a hole in the middle to allow normal plant growth, thus reducing 
the evaporation of water from the soil surface. Throughout the 
experiments, the amounts of water in the pots were determined by 
weighing the pots every day. This weighing recorded the amount 
of daily water loss, corresponding to the daily transpiration of the 
plants. For each pot, at the end of the experiment, cumulative 
daily transpiration was called CWT (the cumulative water 
transpired). Biomass was separated into leaves and stems at 
harvest. Total dry aerial biomass (BM) was obtained after drying at 
80°G for 48 h. WUE was determined at the end of the 
experiment, defined as the ratio of BM to CWT. In addition, a 
dual drought scenario strategy for the two experiments (explained 
in detail below) was studied. 

Experiment conducted in 201 1: scenario of progressive 
water stress 

A randomized complete block design with three replicates was 
used for the progressive water stress treatments (three repli- 
cates x 150 genotypes = 450 plants; called WS). There was another 
replicate (150 plants) that was considered as a well-watered 
treatment, called WW. 

At 1 day after emergence (DAE), 1 7 days after sowing (DAS), all 
600 pots were watered to field capacity, corresponding to 39.5% 
of soil water content (SWC). These 600 pots (WW and WS) were 
kept without irrigation until 17 DAE (Fig. 1A). In these conditions, 
stomatal conductance of the plant was still not affected. We 
calculated that stomatal conductance started to decrease at an 
average SWC of about 21% (unpublished data). 

Starting at 17 DAE, when genotypes reached around 23% of 
SWC, we irrigated the WW treatment to 30% of SWC and we 
maintained this SWC by daily irrigation. The WS treatment was 
kept without irrigation until harvest (during 15 days). 

Two determinations of WUE were made. The first was the total 
water use efficiency, WUE T2 on, calculated by dividing the BM by 
the CWT 31d . CWT 31d is the cumulative water transpired during 
31 days (from 1 to 31 DAE). The second calculation of WUE was 
made during the period when the two treatments differed in their 
soil water content (WW and WS), from 17 to 31 DAE, and called 
WUE E2 on (water use efficiency "estimation"). WUE E2011 was 
calculated by dividing the "estimated biomass" (BM E ), by the 
CWT 15d , calculated from 17 to 31 DAE. BM E = BM - BM 17 , 
where BM 17 is the biomass estimated at 17 DAE. In addition, the 
BM 17 was calculated as follows: BM [7 = (LA 17 /LA :il ) xBM, where 
LA 17 and LA 31 are the leaf areas measured on 17 and 31 DAE, 
respectively. 

Experiment conducted in 2012: scenario of stable SWC 

A randomized complete block design with two treatments and 
two replicates was performed (300 pots per treatment). Treatments 



consisted of two levels of stable SWC which was imposed: well- 
watered (30% of SWC, namely WW) and water-stressed (16% of 
SWC, namely WS) (Fig. IB). 

At 1 DAE (19 DAS), stable water contents corresponding to 
30% of SWC (WW) and 16% of SWC (WS) were maintained for 
23 days (Fig. IB). WUE was calculated by dividing the BM by the 
CWT 23d (WUE T2012 ), where CWT 23d is the cumulative water 
transpired during 23 days (from 1 to 23 DAE). 

Determination of carbon isotope discrimination (CID) 

Carbon isotope composition (8) was calculated relative to the 
international Pee Dee Belemnite (PDB) standard [26] : 8 p i ant = (R^ 
- R sd )/R sd xl000 [%o] where R sa and R sd are the l:5 C: 12 C ratios of 
the sample and the standard, respectively [27]. Carbon isotope 
discrimination (CID), a factor related to isotope fractionation by 
the photosynthetic process relative to the source carbon was then 
estimated as CID = (8 air - 8 p i ant )/(l+ 8 p i ant / 1000) where 8^ is the 
C composition of atmospheric C0 2 , which is assumed to be — 
8.0%o [26]. Before calculating CID, oven-dried leaves of each 
plant were ground into a homogenous fine powder and 2-3 mg 
subsamples were weighed and placed into tin capsules (Elemental 
Microanalysis, UK) to be analyzed using a continuous flow Isotope 
Ratio Mass Spectrometer (Sercon Ltd., Cheshire, UK) at UC 
Davis Stable Isotope Facility (California, USA). 

Genetic map construction 

A set of 9832 SNPs were used to produce an Infinium HD 
iSelect BeadChip (Infinium). These SNPs were selected from 
either genomic re-sequencing or transcriptomic experiments. The 
gDNA from the INEDI RILs population obtained from the cross 
between XRQand PSC8 lines (210 samples) were genotyped with 
the Infinium array. All genotyping experiments were performed by 
Integragen (IntegraGen SA, Genopole Campus 1 - Genavenir 8, 
5 rue Henri Desbrueres, 91000 Evry, France) and the genotypic 
data were obtained with the Genome Studio software (Illumina) 
with automatic and manual calling. From the 9832 SNPs, 7094 
were technically functional with more than 200 samples having a 
genotyping data. From this set of 7094 markers, 2576 were 
polymorphic between XRQ and PSC8 and 2164 did not show 
distortion of segregation in the population. We used CarthaGene 
vl.3 [34] to build the genetic maps. We added the genotypic data 
of markers from a consensus map [35] to the set of the 2164 SNPs 
to assign them to the appropriate LG to the group 0.3 8 in 
CarthaGene. They were ordered using the Ikh 1 -1 function in 
CarthaGene for each group. The genetic map consisted of 2610 
markers located on the 17 LG for a total genetic distance of 
1863.1 cM and grouped on 999 different loci. All data will be 
available through the www.heliagene.org portal. 

Statistical and QTL analysis 

The data were first tested for normal distribution with the 
Kolmogorov-Smirnov test. These data were subjected to analysis 
of variance (ANOVA) and phenotypic correlation analysis 
(Pearson's correlation) using the software of statistical package 
PASW statistics 18 (IBM, New York, USA). Means were 
compared using a Student-Newman-Keuls (SNK) test (P<0.05). 
The broad sense heritability {h ) was then computed from the 
estimates of genetic (cT 2 g) and residual (<r 2 e) variances derived from 
the expected mean squares of the analyses of variance as h 2 = o 2 g/ 
(a 2 g+ a 2 e/ r), where r was the number of replicates. 

QTL identification was performed using MCQTL, software for 
QTL analysis (http://carlit.toulouse.inra.fr/MCQTL/). The 
MCQTL software package can be used to perform QTL mapping 
in a multi-cross design. It allows the analysis of the usual 
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populations derived from inbred lines [28]. MCQTL package is 
comprised of three software applications. The first component, 
TranslateData reads data from MAPMAKER [29] like files. The 
second component, ProbaPop computes QTL genotype probabil- 
ities given marker information at each chromosome location for 
each family and stores them in XML formatted files. The last 
component, Multipop builds the pooled model using the genotype 
probabilities, computes Fisher tests and estimates the model 
parameters [28]. The statistical significance of QTLs was assessed 
using the MCQTL test, which is equal to -log(P-value (F-test)), as 
described in the MCQTL user guide. 

Significant thresholds (P<0.05) for QTL detection were 
calculated for each dataset using 1000 permutations [30] and a 
genome-wide error rate of 0.01 (Type I error). The corresponding 
type I error rate at the whole-genome level was calculated as a 
function of the overall number of markers in the map and the 
number of markers in each linkage group [31]. In our analysis, the 
threshold for the Fisher test (-log(P-value (F-test))) was 3.69 for 
both experiments. This threshold was an average of several 
thresholds of the traits at a significance level of 5% and was 
determined after 1000 permutations. 

In each experiment, the QTL detection was also performed to 
identify QTL for the phenotypic response (called "response 
QTL"), calculated as the difference between two different water 
treatments (WW and WS). This allowed us to detect chromosome 
regions having quantitative effects on traits, depending on the 
environment [32-33]. 

Results 

Genotypic variability and phenotypic correlation 
between water use efficiency (WUE) and carbon isotope 
discrimination (CID) 

In general, a normal distribution was observed for WUE and 
CID traits across the two experiments and water treatments, 
except for WUE T2 oi2 and CID in Exp. 2012 at WW conditions, 
the distributions deviate from normality according to the 
Kolmogorov-Smirnov test (Fig. 2 and 3). As normalizing data 
through transformation may misrepresent differences among 
individuals by pulling skewed tails toward the center of the 
distribution [30], all phenotypic analyses were performed on 
untransformed data. 

Higher mean values for WUE for WS (2.31 to 3.06 g.kg _1 )than 
for WW (1.91 to 2.95 g.kg" 1 ) (Table SI and S2) were observed in 
each experiment. In contrast, higher mean values for CID for WW 
than for WS were also observed in each experiment. In addition, a 
similar range of WUE and CID values was observed in both 
experiments for both WW and WS (for WUE in Exp. 201 1 was 
represented by the WUE E2 on). In addition, significant genotypic 
effects were detected for all traits in Exp. 2011 (Table SI), and 
significant genotypic and SWC effects were detected for all traits in 
Exp. 2012 (Table 1). 

The heritabilities of CID were usually higher than those of 
WUE in both experiments (CID with WUE T2011 or WUE T2012 ), 
except that the heritability of WUE E2 on was higher than that of 
CID (Table SI and 1). 

Significant negative correlations were observed between WUE 
and CID in both experiments (r p = -0.197, P<0.05; r p = -0.409, 
P<0.001; r p = -0.565, P<0.001 for the correlations ofWUE T2011 , 
WUE E2011 , WUE X2012 with the CID, respectively; Fig. 4, Table 
S3, S4 and S5). However, when we determined the correlation 
between WUE and CID for each treatment, we observed a 
positive correlation between the WUE T2 on and CID in Exp. 201 1 
for WS (Fig. 4 and Table S4). In addition, a significant phenotypic 



correlation was observed between Exp. 2011 and 2012 for both 
WUE and CID (Fig. 5). 

QTL identified for water use efficiency (WUE) 

In Exp. 201 1, two QTL for WUE T2011 were detected for WW 
and four QTL for WUE E2011 were detected for WS (Table 2). For 
WW, the QTL were located on LG06 and LG1 1 with the highest 
likelihood odds ratio (LOD) value at 3 cM (QTL of 
WUEllww.11.1) (Fig. SI). The marker for the QTL of 
WUEUww.11.1 was identified between the markers of 
HA005673_395 and HA006174_145 (Fig. 6). For WS, the QTL 
were located on chromosomes LG03 and LG16 (two QTL for 
each chromosome) with the highest LOD value at 6 cM, the QTL 
of WUEellws.16.2, and the marker of this QTL was 
HA017124_226. A "response QTL" for WUE (WUE1 l£ff.06.2) 
was collocated with QTL of WUEllww.06.1. In addition, two 
other "response QTL" were found on LG05 and LG06. The 
additive effects of the WUEllww.06.1 and WUEllww.11.1 were 
-0.14 and 0.11 while the additive effects of the WUEel lws.03.1, 
WUEel lws.03.2, WUEel lws. 16.1 , and WUEellws.16.2 were - 
0.13, 0.13, 0.38 and -0.44, respectively. 

In Exp. 2012, two QTL for WUE T2012 were detected at WW 
and one QTL for WUE T2012 at WS (Table 3). For WW, the QTL 
were detected on chromosome LG13 and LG15 with the highest 
LOD value at 25 cM, the QTL of WUE12ww.l3.1, and the 
markers for this QTL was restor (Fig. 6, Fig. SI). For WS, a QTL 
was detected on chromosome LG09 (QTL of WUE12ws.09.1) with 
the LOD value at 3 cM. The marker for the QTL of 
WUE12ws.()9.1 was identified between the markers of SSL053 
and HA013641_506. In addition, a "response QTL" for WUE 
{WUE12diff.l3.1) was co-located with the QTL of WUE12ww.l3.1 
and CID12ww.l3.1 (Table 4). The additive effects of 
WUE12ws.09.1, WUE12ww.l3.1 and WUE12ww.l5.1 were 0.20, 
0.04 and —0.06, respectively. 

QTL identified for carbon isotope discrimination (CID) 

In Exp. 2011, two QTL for CID were detected at WW and 
three QTL for CID were detected at WS (Table 2). For WW, the 
QTL were located on the same chromosomes of LG06 with the 
highest LOD value at 4.5 cM, QTL of CID1 lww.06.1, and the 
marker of this QTL was ORS483 (Fig. 6, Fig. S2). For WS, the 
QTL were identified on chromosomes LG03, LG06 and LG13 
with the highest LOD value at 5.5 cM, the QTL of CID1 lws.03.1, 
and the marker of this QTL was HA013974_334. Besides, there 
was one "response QTL" detected for CID on chromosome LG02 
(CIDlldiff.02.1) (Table 4). The additive effects were -0.15 and 
0.12 (for QTL of CID1 lww.06.1 and CID1 lww.06.2) while the 
additive effects were -0.13, -0.10, -0.13 (for the QTL of 
Cim lws.03.1, CIDllws.06.1 and CID1 lws. 13.1) (Table 2). 

In Exp. 2012, two QTL for CID were detected at WW and one 
QTL for CID at WS (Table 3). For WW, the QTL were found on 
chromosomes LG13 and LG15 with the highest LOD value of 
8.5 cM, the QTL of CID12ww.l3.1, and the marker for this QTL 
was restor (Fig. 6, Fig. S2). For WS, a QTL was found on 
chromosome LG13 with an LOD value of 2.5 cM; the QTL of 
CID12ws.l3.1, and the marker for this QTL was HACG0018_- 
Contig_l_130. The additive effects for CID12ww.l3.1 and 
CID12ww.l5.1 were 0.20 and 0.07, respectively. The additive 
effect of the QTL of CID at WS {CID12ws.l3.1) was 0.14. 
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Figure 2. Frequency distribution for water use efficiency (WUE) in Exp. 2011 and 2012 of 150 recombinant inbred lines (RILs). 

WUE T2 oii: total water use efficiency "total" in Exp. 2011; WUE E20qq : water use efficiency "estimation" in Exp. 2011; WUE T2 oi2 : water use efficiency 
"total" in Exp. 2012. WW: well-watered; WS: water-stressed. For WUE T2011 and WUE E2011 at WW, data represent 150 RILs (n = 150); for WUE T2011 and 
WUE E20 ii at WS, data represent mean of three replicates of 150 RILs (n = 150); for WUE T20 i 2 at WW and WS, data represent mean of two replicates of 
150 RILs (n = 150). SD: standard deviation. 
doi:10.1371/journal.pone.0101218.g002 
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Table 1. Heritability (/i 2 ) and mean square (MS) of analysis of variance (ANOVA) for water use efficiency (WUE), carbon isotope 
discrimination (CID), biomass (BM) and cumulative water transpired (CWT) for 150 recombinant inbred lines (RILs), two stable soil 
water contents (SWC) and two replicates in Exp. 2012 (n = 600). 





Trait 


h 2 


MS 






Genotype 


Soil water content 


Genotype xsoil water content 


WUE T2012 


0.26 


0.50*** 


28*** 


0.25 ns 


CID 


0.41 


1 .68*** 


1100*** 


0.53 ns 


BM 


0.36 


0.51*** 


180*** 


0.29** 


CWT 23d 


0.36 


40862*** 


31746440*** 


25565*** 



"Significant at P-C0.01, 
♦••significant at P<0.001. 
ns Not significant. 

doi:1 0.1 371 /journal.pone.01 01 21 8.t001 
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Figure 4. Relationship between water use efficiency (WUE) and carbon isotope discrimination (CID) of 150 recombinant inbred 
lines (RILs) in Exp. 2011 and Exp. 2012. Relationship between (A) WUE T20ln and CID in Exp. 2011, (B) WUE T20nl and CID at WW in Exp. 2011, (C) 
WUE T2011 and CID at WS in Exp. 201 1, (D) WUE E2011 and CID in Exp. 201 1, (E) WUE E2011 and CID at WW in Exp. 201 1, (F) WUE E2011 and CID at WS in Exp. 
201 1, (G) WUE T20 i2 and CID in Exp. 2012, (H) WUE T20 i2 and CID at WW in Exp. 2012; (I) WUE T20 i2 and CID at WS in Exp. 2012. Phenotypic correlation (r p ) 
value is provided in each graph. 
doi:1 0.1 371 /journal.pone.01 01 21 8.g004 
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Figure 5. Relationship between (A, B) WUE and (B) CID values for 150 recombinant inbred lines (RILs) determined in two separate 
experiments (Exp. 2011 and 2012). For each trait and experiment, mean of well-watered (WW) and water-stressed (WS) plants were grouped 
together (n = 300). Phenotypic correlation (r p ) value is provided in each graph. 
doi:10.1371/journal.pone.0101218.g005 



QTL identified for related traits: biomass (BM) and 
cumulative water transpired (CWT) 

In Exp. 201 1, three significant QTL for BM, and one QTL for 
each of BM E and CWT 31d at WS were identified (Table 2). These 
QTL were detected on chromosomes LG14, LG15, LG17, LG01 
and LG11. There were only two "response QTL" detected for 
each of BM E and CWT 15d . These QTL were detected on the same 
chromosome, LG06. 

In Exp. 2012, seven QTL were identified for BM under both 
levels of SWC For CWT 2 3 d , five significant QTL were detected 
under both levels of SWC. Further, six "response QTL" for BM 
and CWT 2 3d were identified on chromosomes LG06, LG09, 
LG13 and LG15. 

Discussion 

Genetic variation and relationship between WUE and CID 

In our experiments, increasing drought lead to an increase in 
WUE and a decrease in CID. This result was previously reported 
by Lauteri et al. [36] in sunflower and is well known in other 
crops, such as durum wheat (Triticum turgidum L.) [37], rice (Oryza 
sativa L.) [38] and eucalyptus (Eucalyptus microtheca) [39]. In 
addition, a similar range of values for WUE and CID was 
observed in the two experiments even though their water stress 
patterns differed. That was likely because the population had been 
constructed from parents that had specific responses in non-limited 
and limited water availability [40-41]. From the phenotypic data, 



XRQ exhibited low WUE while PSC8 exhibited high WUE 
(unpublished data). 

CID is highly heritable trait and its heritability is usually higher 
rather than WUE [7,11]. Nevertheless, in the present study, both 
of CID and WUE were influenced by environmental variation 
because the heritability values were below 50% [24]. A previous 
study [42] has shown that heritabilities for CID, measured on 
detached sunflower leaves, were above 50% (74-96%), indicating 
that genetic variance for CID was dominant. However, this result 
was obtained for plants grown in optimal watering conditions. 
Consequendy, CID appeared dependent on genetic and environ- 
mental control. This trait is genetically complex [43], and its 
expression in leaves and other plant tissues varies with the water 
supply. In drought conditions, Rebetzke et al. [24] reported that 
low soil water availability decreases stomatal conductance, which 
can reduce genetic variance and heritability of CID. 

Our work demonstrated the clear relationship between WUE 
and CID in different water regimes. For each water regime and all 
genotypes, we observed negative correlations between WUE and 
CID. These results are in accordance with those of previous work 
in sunflower [36,42], and with those of numerous authors working 
on other crops [6,8,44-47]. In one case of progressive water stress, 
WUE T2 on and CID, were positively correlated. This was 
probably due to the high variability of the soil water content 
during the progressive drought establishment (SWC was gradually 
decreased). A similar result was reported on alfalfa genotypes [48]: 
WUE (mg of dry matter per g H 2 0) was positively correlated with 
CID for plants subjected to progressive water stress during 7 days. 
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Figure 6. Genetic locations of QTL for water use efficiency (WUE), carbon isotope discrimination (CID), biomass (BM) and 
cumulative water transpired (CWT) in the progressive stress experiment (201 1 ) and the stable stress experiment (201 2). Numbers on 
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the left of linkage groups (LG) indicate the cumulative distance in centimorgan (cM) to the first marker at the top LG. Marker names and QTL are 
specified to the right of LG. The same QTLs which are found in a LG are shown in bold. Not all these chromosomes contain the complete markers 
(each chromosome has only been provided by the markers at the top, middle and bottom of LG as well as the markers for identified QTLs). QTL 
confidence intervals were estimated using the two-LOD confidence region. 
doi:10.1371/journal.pone.0101218.g006 



In the WW treatment, the high WUE was correlated with high 
BM and high CWT, while for the WS treatment the high WUE 
was still correlated to high BM but with low CWT. If increase in 
WUE is associated with reduced transpiration, such genotypes are 
often referred to as "conductance type". On the other hand, if 
increase in WUE is correlated with increased photosynthesis, such 
genotypes can be categorized as "capacity types" [49-50]. 
Accordingly, the sunflower genotypes in our study can be 
categorized as an intermediate between "conductance" and 
"capacity" type, unlike rice genotypes that have been categorized 
as "conductance type" [51]. In addition, our results were in 
agreement with several authors [39,52-53] who have suggested 
that plants that use water more efficiently by producing greater 
biomass for a given quantity of water transpired would grow more 
rapidly, resulting in a positive correlation between WUE and 
biomass production. 

QTL identified for WUE and CID 

Our study is the first to identify QTL for WUE and CID in 
sunflower subjected to drought. In Exp. 2011, significant regions 
affecting WUE were identified on four different chromosomes 
(LG03, LG06, LG11, LG16) in two water treatments and 
significant regions affecting CID were identified on three different 
chromosomes (LG03, LG06, LG13) for the same two water 
treatments. From these QTL, we observed a decrease and an 
increase of additive effects (XRQ, indicating that genes having 
both negative and positive effects had been involved in the 
difference in WUE and CID between the parental lines [54]. In 
Exp. 2012, the QTL for WUE were detected on three different 
chromosomes in two water treatments (LG09, LG13 and LG15) 
and the QTL for CID were identified on two different 
chromosomes in these two water treatments (LG13 and LG15). 
All these QTL increased the values of additive effects except the 
QTL of WUE12ws.09.1, indicating that XRQ allele increased the 
traits. These findings provide an explanation for the underlying 
genetic basis of the transgressive variation observed in the 
segregating population. This is in accordance with the argument 
proposed by Chapman et al. [55] and Vargas et al. [56], namely 
that a given QTL can have positive or negative additive effects, or 
none at all, depending on the drought scenario. 

The WUE and CID were controlled by several QTL with small 
genetic additive (XRQ effects, indicating that WUE and CID 
were genetically complex traits [2,57]. Reports evaluating genetic 
analysis for CID in other crops like soybean [58], cotton [59] and 
rice [54] have identified multiple QTL of smaller effect associated 
with the trait. However, in the present study, the QTL for WUE 
and CID explained 42% and 21% of the highest phenotypic 
variance (R 2 ). These R 2 values are higher than those found by 
previous authors for other crops, for example, rice [10,54], wheat 
[24] and barley [60-61]. 

Expression of QTL for WUE and CID across experiments 
and water treatments 

The locations of QTL might be affected by growth stage [54] 
and/or environmental change [62-63]. In our results, the QTL 
for WUEx2oi i and WUEe 2 oi i were found on chromosomes LG03, 
LG06, LG1 1 and LG16 (under WW and WS), whereas the QTL 
for WUE T2 oi2 were found on chromosomes LG09, LG13 and 



LG15 (under WW and WS). These results showed that the 
expression of QTL for WUE differs with micro-environmental 
variations. This variation can be explained by the different water 
regimes in Exp. 2011 and Exp. 2012. 

When the same mapping population is phenotyped in different 
environments, some QTL could be detected in one environment 
but not in others [63] . Collins et al. [64] noted that QTL can be 
categorized according to the stability of their effects across 
environmental conditions. A "constitutive" QTL is consistendy 
detected across most environments, while an "adaptive" QTL is 
detected only in specific environmental conditions or increases in 
expression with the level of an environmental factor. 

The QTL for CID in Exp. 201 1 were detected on chromosomes 
LG03, LG06 and LG13 (WW and WS), whereas the QTL for 
CID in Exp. 2012 were detected on chromosomes LG13 and 
LG15 (under WW and WS). These results indicate that the 
expression of QTL for CID differs in the two experiments and 
different water regimes. Despite CID variation is influenced by 
stomatal conductance and photosynthetic capacity variations 
[7,37], several QTL of the different water regimes have been 
detected on the same chromosome [65-67]. This was the case in 
our study, where the three QTL for CID of the three different 
water regimes were detected on the same chromosome (LG13). 
Therefore, the QTL for CID in this study can be considered as a 
"constitutive" QTL. Additionally, the constitutive QTL for CID 
was consistent with the result of phenotypic correlation that 
genotypic ranking for this trait was consistently maintained in the 
two experiments. 

Some QTL for WUE and CID and related traits were located 
on the same chromosome or on a similar QTL position (co- 
localization). The QTL for WUE T2 oi 2 for WW {WUE12ww.l3.1 
and WUE12ww.l5.1) had a similar QTL position (26.20 and 
77.10 cM) as the QTL for CID for WW (CID12ww.l3.1, 
CID12ww.l5.1). The QTL for CID (CID12ws.l3.1) for WS was 
associated with the QTL for WUE T2012 for WW 
[WUE 12ww. 13.1). This QTL was detected on chromosome 
LG13 (QTL position: 30.80 cM) near the QTL of CID12ww.l3.1 . 
The occurrence of QTL associated with different traits at the same 
locus may be explained by the fact that (i) the QTL are closely 
linked genetically or (ii) a single locus controls multiple traits and a 
gene may have pleiotropic effects [54] . 

We have observed a common genetic basis for WUE and CID 
in each experiment. Using the same mapping population under 
different water stress treatments helped us to characterize 
consistent genomic region (by QTL). Kiani et al. [68] indicated 
that QTL which was induced only by drought might be associated 
with mechanism(s) of sunflower drought response and they 
proposed that the QTL which can reduce trait difference between 
well-watered and water-stressed conditions should have an effect 
on drought tolerance because of their contribution to trait stability. 
Our study in Exp. 2011 showed that the QTL for CID on 
chromosome LG06 were repeatable across two different water 
treatments (WW and WS). In Exp. 2012, the QTL for CID on 
chromosome LG13 have been repeatable across two different 
water treatments (WW and WS). 

All these QTL which are common across different water 
treatments might be useful for marker-assisted selection (MAS). 
Identification of QTL influencing several traits could increase the 
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efficiency of marker-assisted selection (MAS) and hasten genetic 
progress [69]. Ribaut et al. [70] noticed that in the design of the 
best-possible breeding strategy using MAS, additional traits and 
criteria have to be considered. For each trait of interest, some of 
the criteria are the number of QTL detected, the percentage of 
phenotypic variance that they explain, the total percentage of the 
genome that they represent, and their stability across different 
environments. Regarding these arguments, our study has shown 
that CID is the most interesting trait and should be useful for 
MAS, where three QTL overlapped on chromosome LG06 (CID 
for WW and WS in Exp. 2011), and three QTL across three 
different water treatments were co-localized on chromosome 
LG 1 3 with phenotypic variance (R 2 ) ranges from 7 to 21%. 
Further, these QTL and other co-localized QTL on chromosomes 
LG06 and LG13 were identified in the near-centromeric region 
(inferior to superior position explained from 0 to 60.06 cM, and 
from 0 to 62.45 cM for LG06 and LG13, respectively), because 
those chromosomes are classified as a metacentric type [71-72]. 

Co-localization of QTL for WUE and CID with related traits 

In this study, we also detected QTL for the related traits BM 
and CWT on the same chromosome of the QTL for WUE and/ or 
CID (for WW and WS). These were observed in Exp. 2012, where 
two of four QTL for BM for WW (BM12ww.l3.1 and 
BM12ww.l5.1) were detected on chromosomes LG13 and LG15, 
and co-located with the QTL for WUE T2012 and CID for WW 
(WUE12ww.l3.1, CID12ww.l5.1). For WS, the identifications of 
the QTL for the related traits showed a similar trend. The QTL of 
BM12w.s.l3.1 (QTL position: 21.20 cM) was detected on chro- 
mosome LG13, as the QTL of CIDllws.13.1, CID12ww.l3.1, 
CID12w.s.l3.1 and WUE12ww.l3.1 have been identified. These 
indicated the possibility of genetic association of WUE and CID 
with the accumulation of biomass. Consistent with this, Kiani 
et al. [68] identified a QTL for total dry matter in water-stressed 
conditions on chromosome LG13 using another population of 
sunflower. Interestingly, this QTL overlapped with osmotic 
adjustment, grain yield, and plant height. Thereby the common 
genetic basis for WUE, CID, productivity and osmotic adjustment 
will lead to an improved understanding of drought tolerance 
genes. In addition, evidence of overlapping QTL of productivity 
and osmotic adjustment have been observed by several authors 
[73-75]. However, further study is obviously required to 
determine the genetic control of osmotic adjustment or hydraulic 
conductance and their inter-relationships with WUE and CID. 

For CWT, the QTL of CWT12ww.l5.1 was detected on 
chromosome LG15 with the QTL position at 79.10 cM near the 
marker at position of 77.10 cM where the QTL of 
WUE12ww.l5.1 and CID12ww.l5.1 have been identified. Not far 
from these positions, a QTL of CWT12ws.l5.1 was also detected 
(QTL position: 49.5 cM). These indicated out that the cumulative 
water transpired in WW and WS is genetically and closely related 
with WUE and CID in non-limited water availability. In addition, 
the maintenance of biomass accumulation under stable water 
stress should be considered as an efficiency process between 
transpiration, biomass accumulation and its partitioning between 
non-drought and drought conditions [64] . Therefore, the increase 
in WUE (i.e. the amount of biomass produced per unit of 
transpired water) might seem to be ideal candidate mechanism for 
drought-prone environments. 

Identifying the "response QTL" for WUE and CID 

In our work, we calculated the "response QTL" to provide new 
insight into the genetic architecture of WUE and CID, which, 
unlike a "common" phenotypic trait, is rarely considered in QTL 



analysis. Water use traits and their response are of primary 
importance to plant growth and survival. Although we have a 
growing understanding of the genetic and molecular drivers of 
water use traits and WUE as well as CID, response QTL of those 
traits has received relatively little attention. 

We detected three QTL of "response QTL" for WUE on 
chromosomes LG06 and LG13. From these two chromosomes we 
have also identified the QTL for WUE T2011 and WUE T2 oi2 for 
WW, indicating, at least under the conditions imposed in these 
experiments, that response QTL was controlled by loci that 
determine the main trait value under a specific treatment. This 
was in agreement with Kliebenstein et al. [76-77] who evaluated 
the response QTL between control and methyl jasmonate (Meja)- 
treated plants of Arabidopsis thaliana. They reported that significant 
QTL that influenced response between control and Meja-treated 
plants also affected the main trait value in at least one of the two 
environments, which was called the "allelic sensitivity" model. 

In contrast, an independent response QTL, was also observed 
for several traits, for example the response QTL for WUE T2011 on 
chromosome LG15 {WUE1 ldiff.05.1), CID on chromosome LG02 
(CIDdiffl 1.02.1), and CWT 23d on chromosome LG06 
(CWT12diff.06.1). This observation was not consistent with 
Kliebenstein et al. [77], however, it was in agreement with an 
argument of Schlichting and Pigliucci [78] who suggested the 
"gene regulation" model must exist, and is not always controlled 
by loci that are expressed within at least one of the two 
environments. 

As for the prospects for these aspects, characterization of the 
genes underlying QTL that control the differential WUE and CID 
regulation might generate a detailed understanding of the 
molecular and biochemical basis for water use traits in sunflower 
and how this alters phenotypic response in more complex 
environments. 

Importance of high WUE or low CID for sunflower 
breeding: use of the identified markers for MAS 

This is the first genetic quantitative analysis and QTL mapping 
for WUE and CID in sunflower. We investigated two drought 
scenarios and evaluated genetic variation of sunflower lines to 
identify genetic control and physiological processes that could 
explain genotypic differences in the response to drought stress. 
The present study proved that, in sunflower, selection for CID can 
be considered in initial screening to improve WUE. However, this 
merits further investigation in other populations. 

Many QTL (particularly for CID) have been reported in the 
literature. However, very few with large effects have been 
adequately exploited in crop breeding programs. The majority 
of the favorable alleles for identified QTL are to be found in 
journals on library shelves rather than in crop cultivars improved 
by introgression or selection of these favorable QTL alleles [79]. 
Nevertheless, Condon et al. [80] reported the release of a new 
high-yielding wheat variety in droughted environments after a 
breeding process in which selection for low CID in non-droughted 
plants led to high WUE. 

In conclusion, our results emphasize that the near-centromeric 
region of chromosomes LG06 and LG13 are a "reliable" region 
for MAS due to the co-localization of the QTL for CID with 
several QTL for WUE, BM and CWT. Indeed, the best strategy 
for using molecular markers should combine selection for QTL 
involved in the expression of CID. 

This paper complements the study of Vincourt et al. [25] and 
Rengel et al. [81] that exploited the INEDI RIL population in 
analyzing genetic variation of agronomic and physiological traits, 
making it possible to establish strategies for a sunflower breeding 
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program and provide a basis for identification of the molecular 
components of a genotype x environment interaction. 
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